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We present measurements of the differential cross section dcr/dp^ for the inclusive production of 
a photon in association with a fe-quark jet for photons with rapidities \y''\ < 1.0 and 30 < p^ < 
300 GeV, as well as for photons with 1.5 < \y''\ < 2.5 and 30 < pj < 200 GeV, where p^ is the 
photon transverse momentum. The 6-quark jets are required to have pr > 15 GeV and rapidity 



Jet I 



< 1.5. The results are based on data corresponding to an integrated luminosity of 8.7 fb 



recorded with the DO detector at the Fermilab Tevatron pp Collider at ^/s =1.96 TeV. The measured 
cross sections are compared with next-to-leading order perturbative QCD calculations using different 
sets of parton distribution functions as well as to predictions based on the fcx-factorization QCD 
approach, and those from the SHERPA and pythia Monte Carlo event generators. 



PACS numbers: 13.85.Qk, 12.38.Bx, 12.38.Qk 



In hadron-hadron collisions, high-energy photons (7) 
emerge unaltered from the hard scattering process of two 
partons and therefore provide a clean probe of the parton 
level dynamics. Study of such photons (called prompt or 
direct) produced in association with a 6-quark also pro- 
vides information about the 5-quark and gluon (g) parton 
distribution functions (PDFs) of the incoming hadrons. 
Such events are produced in Quantum Chromodynam- 
ics (QCD) primarily through the Compton-likc scattering 
process gb — >■ 76, which dominates up to photon trans- 
verse momenta (pj) of sa 70 GeV, and through quark- 
antiquark annihilation (7^ — >■ 75 — >■ 766, which dominates 
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at high p^ m . The inclusive 7 + 6 production may also 
originate from partonic processes like gg — > bb or bg — > bg, 
where the final state 6-quark or gluon fragments into 
a photon [l|. However, photon isolation requirements 
substantially reduce the contributions from this process. 
The measurements of the differential cross section as a 
function of p^ and the photon (and/or 6-jet) rapidity 
can be used to test the 7 + 6 production mechanism and 
the underlying dynamics of QCD hard-scattering subpro- 
cesses with different momentum transfer scales Q and 
parton momentum fraction x. Measurements involving 
7/Z-boson and 6-quark final states have previouslybeen 
performed by the DO and CDF collaborations [21461. In 
comparison to the previous 7 + 6 measurement [6|, we 
now consider not only the leading (in pt) 6-jet, but all 
6-jets in the event. To increase statistics in p^ bins, we 
have also extended the |y''°*| region which results in a 
larger contribution from the annihilation process. The 
large integrated luminosity recorded with the DO detec- 
tor in pp collisions at ^/s = 1.96 TeV at the Fermilab 



Tevatron Collider and more advanced photon and 6-jet 
identification tools [THSl enable us to perform more pre- 
cise measurements and to extend them in kinematic re- 
gions previously unexplored. 

In this Letter, we present measurements of the inclu- 
sive 7 -I- 6-jet production cross sections using data col- 
lected from June 2006 to September 2011. The cross 
sections are measured as a function of p^ in the photon 
rapidity regions, \y^\ < 1.0 (central) and 1.5 < \y^\ < 2.5 
(forward). The rapidity, y, is related to the polar scat- 
tering angle 6 with respect to the proton beam axis by 
y = ^ ln[(l -t- Pcos6)/{l ~ f3cos9)], where /3 is defined 
as the ratio between momentum and energy /3 ~ \p\/E. 
The photons are required to have 30 < p^ < 300 GcV 
in the central rapidity region and 30 < pj < 200 GeV 
in the forward region. The 6-jets are required to be 
within Iw-^'^'l < 1.5 and to have transverse momentum 



p!^ > 15 GeV. This allows us to probe the dynamics 
of the production process in a wide kinematic range, not 
studied before in other measurements of a vector boson 
-|- 6-jet final state. The measurement covers parton mo- 
mentum fractions in the range 0.007 < a; < 0.4. Figure 
[T] shows the fractional contributions of the gb — s- 76 sub- 
processes to the total cross section of 7 -f- & production 
with photons in the central and forward photon rapid- 
ity regions as a function oi p^. The curves are obtained 
using signal processes gb — > 76 and qq — > 766 simulated 
with the PYTHIA event generator [lO| . It can be seen that 
the Compton-like contribution is large at small p^ and 
decreases with growing p^, with the annihilation process 
contribution having the opposite behavior. 
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FIG. 1: Fractional contribution of the gb — )> 76 subprocess 
to the associated production of direct photon and fe-jet as 
a function of pj in the events with photons in the central 
and forward rapidity regions. The fractions are calculated 
using PYTHIA 6.4 [10) and the CTEQ6.1L parton distribution 
functions [Till. 



The DO detector is a general purpose detector dis- 
cussed in detail elsewhere Il21. The subdctectors most 



relevant to this analysis are the central tracking sys- 
tem, composed of a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT) embedded in a 1.9 T 
solenoidal magnetic field, the central preshower detector 
(CPS), and the calorimeter. The GPS is located immedi- 
ately before the inner layer of the central calorimeter and 
is formed of approximately one radiation length of lead 
absorber followed by three layers of scintillating strips. 
The calorimeter consists of a central section with cover- 
age in pseudorapidity of |7?det| < 1-1 [l3l, and two end 
calorimeters covering up to |??det| ~ 4.2. The electro- 
magnetic (EM) section of the calorimeter is segmented 
longitudinally into four layers (EMi, i = 1 — 4), with 
transverse segmentation into cells of size Aj^dot x A0dot = 
0.1 X 0.1 [U, except EMS (near the EM shower maxi- 
mum), where it is 0.05 x 0.05. The calorimeter allows 
for a precise measurement of the energy and direction of 
electrons and photons, providing an energy resolution of 
approximately 4% (3%) at an energy of 30 (100) GcV, 
and an angular resolution of 0.01 radians. The energy 
response of the calorimeter to photons is calibrated us- 
ing electrons from Z boson decays. Since electrons and 
photons shower differently in matter, additional energy 
corrections as a function of y'' are derived using a detailed 
GEANT-based [IJ] simulation of the DO detector response. 
These corrections are largest, « 2%, at photon energies 
of about 30 GeV. The data used in this analysis satisfy 
DO data quality requirements and are collected using a 
combination of triggers requiring a cluster of energy in 
the electromagnetic (EM) calorimeter with loose shower 
shape requirements, and cor resp ond to an integrated lu- 
minosity of 8.7 ±0.5 fb~^ [15|. The trigger efficiency 
is w 96% for photon candidates with p^ ^ 30 GeV and 
« 100% for pI > 40 GeV. 

Offline event selection requires a reconstructed pp in- 
teraction vertex [16[ within 60 cm of the center of the 
detector along the beam axis. The efficiency of the ver- 
tex requirement is « (96 — 98)%, depending on p"^. The 
missing transverse momentum in the event is required to 
be less than 0.7p^ to suppress background from W ^ ev 
decays. Such a requirement is highly efficient for signal 
events, with an efficiency > 98% even for events with 
semi-leptonic heavy-flavor quark decays. 

To reconstruct photon candidates, projective towers 
of calorimeter cells with large deposits of energy are 
used as seeds to create clusters of energy in the EM 
calorimeter in a cone of radius TZ = 0.4, where TZ = 
^J{/S.riY + (A0)^. Once an EM energy cluster is formed, 
the final energy {E-em) is obtained summing the energies 
of all the calorimeter cells in a smaller cone of 7?. = 0.2. 
Photon candidates are required to have: (i) > 97% of 
their energy in the EM section; (ii) calorimeter isola- 
tion X = [^tot(0.4) - ^em(0.2)]/£;em(0.2) < 0.07, where 
Etot(JZ) [£^em(7?.)] is the total [EM only] energy in a cone 
of radius TZ] (iii) scalar sum of pt less than 1.5 GeV, cal- 
culated from all tracks with px > 0.5 GeV originating 



from the pp primary interaction point in an annulus of 
0.05 <TZ < 0.4 around the EM cluster; and (iv) energy- 
weighted EM shower width consistent with that expected 
for an electromagnetic shower. To suppress electrons 
misidentified as photons, the EM clusters are required 
to be not spatially matched to significant tracker activ- 
ity, either a reconstructed track or, in the central rapidity 
region, a density of hits in the SMT and CFT consistent 
with that of an electron [17|. In the following, this re- 
quirement is referred to as the "track-match veto." 

To further suppress jets misidentified as photons, 
an artificial neural network (7-NN) discriminant is de- 
fined [8[. It relies on differences between photons and 
jets in tracker activity, energy deposits in the calorime- 
ter, and in the CPS for the central photons/jets. This 7- 
NN is trained using PYTHIA [ifll Monte Carlo (MC) sam- 
ples of photon and jet production, which are processed 
through a GEANT-based IJ] simulation of the detector 
geometry and response. In order to accurately model 
the effects of multiple pp interactions and detector noise, 
events from random pp crossings with a similar instan- 
taneous luminosity spectrum as in data are overlaid on 
the MC events. These MC events are then processed us- 
ing the same reconstruction code as for the data. The 
7-NN performance is verified using a data sample con- 
sisting of photons radiated from leptons in Z boson de- 
cays {Z -^ £+£--/, i^e.n) [il. The 7-NN output Onn 
distributions for photons in data and MC are in good 
agreement and exhibit a significant separation from the 
distribution for misidentified jets [71 18|. Photon candi- 
dates are required to have Onn > 0.3, which is « 98% 
efficient for photons. 

We calculate corrections to the observed number of 
candidate events to account for the photon detection 
efficiency and for the acceptance of the selection using 
simulated samples of 7 + fe-jet events. In these samples, 
the photon is required to be isolated at particle level by 
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El < 2.5 GeV, where E^^ is the to- 



tal transverse energy of particles within a cone of radius 
TZ = 0.4 centered on the photon, and E^ is the photon 
transverse energy. Here, the particle level includes all 
stable particles as defined in Ref. [19|. Signal events are 
generated using the SHERPA [20| and PYTHIA event gener- 
ators, processed through a GEANT-based [iJI simulation 
and events reconstruction as described above. The accep- 
tance is driven by the selection requirements in 77dot (ap- 
plied to avoid edge effects in the calorimeter regions used 
for the measurement) and (/)dct m the central rapidity 
region (to avoid periodic calorimeter module boundaries 
[12| that bias the EM cluster energy and position mea- 
surements), photon rapidity y'^ and energy, and bin-to- 
bin migration effects due to the finite energy and angular 
resolution of the EM calorimeter. The acceptance varies 
within (82 — 90)% with a relative systematic uncertainty 
of (2 — 5)%. The EM clusters reconstructed in the accep- 
tance region are required to pass the photon identification 



criteria listed in the previous paragraph. Average correc- 
tion factors to account for differences between data and 
simulations are obtained with the SHERPA events, while 
the difference from the corrections obtained with pythia 
is used as systematic uncertainty. Small differences be- 
tween data and MC in the photon selection efficiencies 
are corrected for with suitable scale factors derived us- 
ing control samples of electrons from Z boson decays, as 
well as photons from the radiative Z boson decays [l8|. 
The total efficiency of the above photon selection criteria 
is (68 — 85)%, depending on the p^ and rapidity region. 
The systematic uncertainties on these values are 3% for 
\y''\ < 1.0 and 7.3% for 1.5 < \y''\ < 2.5 and are mainly 
due to uncertainties caused by the track-match veto, iso- 
lation, and the 7-NN requirements. The contamination 
from Z{-^ e"'"e~)-|-jet and VF(— >■ eiy)+]et events is esti- 
mated from the simulation and is found to be negligible 
(< 1%) for both photon rapidity regions. 

At least one jet with pi^* > 15 GeV and lyj***] < 1.5 
must be present in each selected event. Jets are recon- 
structed using the DO Run II algorithm [2l| with a cone 
radius of 7?. = 0.5. The jet acceptance with respect to 
the pl^ and \y^'^^\ kinematic cuts varies between 88% and 
100% in different photon px bins. The uncertainties on 
the acceptance due to the jet energy scale, jet energy 
resolution, and difference in energy scale correction be- 
tween light flavor and ti-jets vary between 1% and 7%, 
increasing for smaller p^. The jet is required to have 
at least two associated tracks with pT > 0.5 GeV with 
at least one hit in the SMT. The track with the high- 
est pt must have pT > 1.0 GeV. These criteria ensure 
that there is sufficient information to classify the jet as 
a heavy-flavor candidate and have a typical efficiency of 
about 90%. Light jets (caused by light quarks or gluons) 
arc suppressed using a dedicated artificial neural network 
(6-NN) [9| that employs the longer lifetimes of heavy- 
flavor hadrons relative to their lighter counterparts. The 
inputs to the 6-NN combine several characteristic quan- 
tities of the jet and associated tracks to provide a con- 
tinuous output value that tends towards one for 6-jets 
and zero for the light jets. The 6-NN input variables 
providing most of the discrimination are the number of 
reconstructed secondary vertices (SV) in the jet, the in- 
variant mass of charged particle tracks associated with 
the SV (Mgv), the number of tracks used to reconstruct 
the SV, the two-dimensional decay length significance of 
the SV in the plane transverse to the beam, a weighted 
combination of the tracks' transverse impact parameter 
significances, and the probability that the tracks associ- 
ated with the jet originate from the pp interaction vertex. 
The jet is required to have a 6-NN output > 0.3. Depend- 
ing on pj, this selection is (40 — 52)% efficient for 6-jets 
with systematic uncertainties of (6 — 23)% for the j + b 
events with {y'^l < 1.0 and of (7 — 11)% for those with 
1.5 < \y'^\ < 2.5, both increasing as a function of pj- 
Only 0.2 — 0.4% of light jets are misidentified as heavy- 



flavor jets, comprising 7% to 10% of the final sample. 

After all selection requirements, 199,515 (139,710) 
events remain in the data samples with the central (for- 
ward) photons. In addition to events with light- flavor jets 
a second source of background is represented by multi-jet 
events in which one jet is misidentified as a photon. To 
estimate the photon purity, the 7-NN distribution in data 
is fitted to a linear combination of templates for photons 
and jets obtained from simulated 7 -I- jet and dijet sam- 
ples, respectively. An independent fit is performed in 
each p^ bin, yielding photon purities between 62% and 
99% for the events with the central photons and between 
40% and 55% for the events with the forward photons. 
The obtained photon fractions are shown in Fig. [2] The 
p^ dependence of the purity is fitted in each region using 
a two-parameter function V = l — exjp{a + bp^). The sys- 
tematic uncertainties on the fit are estimated using two 
alternative fitting functions. These photon purities differ 
at most by 7% when compared with those obtained for 
inclusive events, i.e. without the requirement of a heavy- 
flavor jet. An additional systematic uncertainty in the 
photon fractions due to the fragmentation model imple- 
mented in PYTHIA is also taken into account [22] ■ This 
uncertainty is estimated by varying the production rate 
of 7r° and rj mesons by ±50% with respect to their central 
values [231 . It is found to be about 6% at p^ ~ 30 GeV, 
2% at pI ~ 50 GeV, and < 1% at pj > 70 GeV. 

The fractions of 6-jets are determined by fitting Afsv 
templates for b, c and light jets to the data. Jets from 
&-quarks tend to have larger values of Mgy, in contrast 
to light jets. For b- and c-jets, the templates are ob- 
tained from simulation, while the light jet template is 
derived from a data sample, enriched in light jets, re- 
ferred to as negatively tagged data (NT data). The NT 
data comprises the jets that have negative values for some 
of the inputs to the 6-NN algorithm (such as negative de- 
cay length and negative impact parameter significance) 
which are caused by detector resolution effects [9| . After 
correcting the NT data for the small contamination from 
heavy-flavor jets, we have verified that the M^v template 
shapes in NT data and light jets in the MC simulation 
agree well. 

The result of a maximum likelihood fit to Msv tem- 
plates, normalized to the number of events in data, is 
shown in Fig. |3] for central photons with 50 < p^ < 
60 GeV, as an example. As shown in Fig.|4l the estimated 
fraction of &-jets grows with pj, from about 35% to about 
42%. The corresponding relative uncertainties range be- 
tween (4-24)%, increasing at higher p'^ and are domi- 
nated by the limited data statistics. The data corrected 
for the photon and jet acceptance, for reconstruction ef- 
ficiencies, for the contribution of background events, and 
for bin-migration effects, are presented at the particle 
level, as defined in Ref. [19|. 

The differential cross sections of 7 4- 6 production are 
extracted in nine (six) bins of pj for central (forward) 



photons, and are listed in Tables U and [III The results 
arc also shown in Fig. [5] as a function of pj. for the two 
photon rapidity intervals. The data points are plotted 
at the value of pj for which the value of the smooth 
function describing the cross section equals the averaged 
cross section in the bin [2J]. The 7 -I- 5-jet simulated 
sample with Sherpa has been used to determine mean p^ 
values. 

The cross sections with the central (forward) photons 
fall by about four (three) orders of magnitude in the 
range 30 < pj < 300 (200) GeV. The statistical un- 
certainty on the results ranges from 2% in the first p^ 
bin to « 11% in the last pj bins, while the total system- 
atic uncertainty varies between 12% and 36%. The main 
source of the uncertainty at low p^ is due to the photon 
purity (up to 8%), the 6-jet fraction fit (6 — 7%), and the 



luminosity (6.1%) [15|. At higher pi^, the uncertainty is 



dominated by the fractions of 6-jets and their selection 
efficiencies. Systematic uncertainties are highly bin-to- 
bin correlated for the first three p^ bins, while the total 
systematic uncertainty is nearly uncorrelated across the 
bins at p^ > 70 GeV. 

Ncxt-to-leading order (NLO) perturbative QCD pre- 
dictions, with the renormalization scale ^r, factoriza- 
tion scale /iF, and fragmentation scale ^/ all set to pZ, 
are also given in Tables U and HH These predictions [l| 
are based on a phase space slicing method used to cal- 
culate the cross section analytically [25|. The uncer- 
tainty from the choice of scale is estimated through a 
simultaneous variation of all three scales by a factor of 
two, i.e. for hr^fj = O.SjOj. and 2pJ. The predictions 
utihze CTEQ6.6M PDFs [ll[ and are corrected for non- 
perturbative effects of parton-to-hadron fragmentation 
and multiple parton interactions. The latter are eval- 
uated using SHERPA and pythia MC samples using their 
default settings [lO|, l2Q|. The overall correction varies 
from about 0.90 at 30 < p^ < 40 GeV to about 0.95 at 
high p^, and an uncertainty of < 2% is assigned to ac- 
count for the difference between the two MC generators. 

The prediction based on the fcr-factorization ap- 
proach 26|,|27| and unintegrated parton distributions [28[ 
are also given in Tables [J and HIl The /cT-factorization 
formalism contains additional contributions to the cross 
sections due to resummation of gluon radiation diagrams 
with k^ above a scale /z^, 0(1 GeV), where kr denotes 
the transverse momentum of the radiated gluon. Apart 
from this resummation, the non- vanishing transverse mo- 
mentum distribution of the colliding partons are taken 
into account. These effects lead to a broadening of the 
photon transverse momentum distribution in this ap- 
proach [26[ . The scale uncertainties on these predictions 
vary from about 31% (-28%) at 30 < p^^ < 40 GeV to 
about 5% (14%) for the central photons and 26% (-13%) 
for the forward photons in the last p^ bin. 

Tables U and [ll] also contain predictions of the pythia 
[lot MC event generator with the CTEQ6.1L PDF set. It 
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FIG. 2: Photon purity as a function of p^ in the selected data sample in the central rapidity region \y''\ < 1.0 (a) and the 
forward rapidity region 1.5 < \y~'\ < 2.5 (b). 



includes only 2 — > 2 matrix elements (ME) with gb — > "fb 
and qq — >■ 75 scatterings (defined at LO) and, with 
g —i' bb splitting in the parton shower (PS). We also pro- 
vide predictions by the SHERPA MC event generator [20[ 
with the CTEQ6.6M PDF set |Il|. For 7 + 6 production, 



SHERPA includes all the MEs with one photon and up to 
three jets, with at least one b-ict in our kinematic region. 
In particular, it accounts for an additional hard jet that 
accompanies the photon associated with a bb pair. Com- 
pared to an NLO calculation, there is an additional bene- 
fit of imposing resummation (further emissions) through 
the consistent combination with the PS. Matching be- 
tween the ME partons and the PS jets follows the pre- 
scription given in Ref. [29|. Systematic uncertainties are 
estimated by varying the ME-PS matching scale by ±5 
GeV around the chosen central value [30[. As a result, 
the SHERPA cross sections vary up to ±7% for the central 
photons and up to — 25%/+17% for the forward photons, 
where the uncertainty is largest in the first p^ bin. 

All the theoretical predictions are obtained including 
the isolation requirement on the photons E!f° < 2.5 GeV 
at the particle level. The predictions arc compared to 
data in Fig. [5] as a function of p^. The ratios of data 
over the NLO QCD calculations and of different QCD 
predictions or MC simulation to the same NLO QCD 
calculations are shown in Fig. [5] as a function of pj for 
both central and forward photons. The ratio of NLO 
predictions with CTEQ6.6M to those with mstw2008 [3i| 
and ABKm09nlo [34I are also shown on the plots. 



The measured cross sections are in agreement with the 
NLO QCD predictions within theoretical and experimen- 
tal uncertainties in the region up to pj < 70 GeV, but 
show notable disagreement for larger p^. The cross sec- 



tion slopes in data significantly differ from the predicted 
ones in both photon rapidity regions. The results indi- 
cate a need for higher order perturbative QCD correc- 
tions in the large p^ region, that is dominated by the 
annihilation process qq -^ "fg {g ^- bb) , and resummation 
of diagrams with additional gluon radiation. The QCD 
predictions from the fcT-factorization approach is in bet- 
ter agreement with data. The best agreement is obtained 
with the SHERPA MC that allows up to two extra hard 
partons (jets) in ME in addition to the 6-quark (jet) and 
includes the consistent treatment of the possible contri- 
butions from the parton shower. However, correcting the 
SHERPA predictions by including additional higher order 
contributions, as is done, for instance, for VF-|-jets events, 
would be desirable (33| . 

In conclusion, we have performed a measurement of the 
differential cross section of inclusive production of photon 
in association with 6-jets at the Tevatron pp collider. The 
results cover the kinematic ranges 30 < pj < 300 GeV 
with |j/T| < 1.0, and 30 < p^^ < 200 GeV with 1.5 < 
\y''\ < 2.5. A good description of the data can only be 
achieved by including higher order corrections into the 
NLO QCD predictions, which are currently present as 
additional real emissions in the SHERPA MC generator. 
These results can be used to improve the description of 
background processes in searches for the Higgs boson or 
for new phenomena beyond the SM at the Tevatron and 
the LHC in final states involving the production of vector 
bosons in association with 6-jets. 

We are grateful to the authors of the theoretical cal- 
culations, T. Stavreva, J. Owens, N. Zotov, S. Schumann 
and F. Siegert, for providing predictions and for many 
useful discussions. 
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TABLE I: The 7 + 6 production cross section da/dpy in bins of p^ for \y~'\ < 1.0 together with statistical, <5(Tstat, and systematic, 
Sasyst, uncertainties. The last four columns show theoretical predictions obtained within NLO QCD, Ict factorization, pythia 
and SHERPA event generators. 
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TABLE IL The 7 + 6 production cross section da/dp^ in bins of p^ for 1.5 < |j/^| < 2.5 together with statistical, Sagtat, 
and systematic, (5o"ayat, uncertainties. The last four columns show theoretical predictions obtained within NLO QCD, ki 
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FIG. 6: The ratio of 7 + 6 production differential cross sections between data and NLO QCD predictions with uncertainties 
for the rapidity regions \y^\ < 1.0 (a) and 1.5 < \y~'\ < 2.5 (b). The uncertainties on the data include both statistical (inner 
error bar) and full uncertainties (entire error bar). Also shown are the uncertainties on the theoretical QCD scales and the 
CTEQ6.6M PDFs. The ratio of NLO predictions with CTEQ6.6M to those with MSTW2008 [sJl and ABKm09nlo [sl] are also 
shown. 



